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ORIGINAL PAGE Td * 


a.l. CARETS — Land use analysis : No change, 

a. 2. CAP.STS — Land use clirriatology : Our first S- 192. computer tapes 

arrived during the period of this report. We have r.adc considerable 
progress in the reading and manipulation of the S-192 type data', \^*ith 
the assistance and advice of Mr. Robert M^^^elly of the USGS Computer 
Center. We first ran print-out naps from channel 21 (thermal band, 
low sampling rate) of the Baltimore-Wa shine ton area where our prime 
test site for the climatology project is located. The first few print' 
outs V7cre block smoothed through a 10 x 10 filter. This facilitated 
the printing of the entire scan width on a single print-out sheet. 

More recently v;e have received an unsmoothed print-out of the first 
one-hundred scan-line pixels along a track form Liberty Reservoir on 
the west to the Chesapeake Bay on the cast. No analysis of this 
scene, however, has taken place, 

A preliminary list of data manipulation and output recuirements for 
S-192 tapes has been compiled. The primary constituents are: 

1) The computation of range, mean standard deviation and 
coefficient of sk.cvniess for all pixelr, within a single 
spectral band and specific area of interest. 

2) Histogram plots of pixel values in PC|1 count intervals 
versus nuirber of occurences in each interval, vitViin a 
single spectral band, within a specific area of interest, 

3) Application of calibration factor, according to a specified 
function, to each data element pixel vrithin a specific 
area of interest. 

A) Production of digital map print-outs t/ith and without block 
smoothing . In the case of smoothed versions a variety of 
matrix operator sixes, e.g. A x A, 5 x 5, 10 x 10, etc., 
would he desirable. 

In the analysis thus far, wq have assumed that the internal scanner 
calibration has already been performed (by }tASA) and that each pixel 
value that we are reading on the tape has already been cprrected as 
a result of that calibration. Ve also assume that no corrections 
have been made for atmospheric water vapor. 

a. 3, Census cities : No change, 

b . RecOTniCiendations concerning decision and/or actions required to 
ensure attainment of the experiment’s scientific objecTiVes i With 

respect to the S055-4 data PCP4 count interval change re<|Ueste<i in the 
last report the nestJ intervals, as recommended, may not be as precise 
as desired. During the coming weeks a more exact determination Of 
these intervals should be made and subsequently forwarded in the 

February monthly report. 



^ • Ex pected accor^;plishH ierLt:s ?, t'lio next report i periods : Conversion 

of rau PCM counts to radiance values and then to tGn])Graturc in ®C 
should be completed by the end of January. VJhcn this is accomplished 
these temperature values uill then be compared 'to the results of 
the Outcalt urban climate simulator, probably here on the USGS computer. 

• Significant results and their r elationship to practical applications 
or operational problems ; 

1. llaps originally created under USGS Contract 14--0o-0001-11914 

at U.C. Pdverside have been refined and a paper lias been prepared 
V7hich describes their making for presentation at the forthcoming 
ASP annual convention in Washington, D.C. These maps, v;hich 
merge, electro-optical scanner technology mith automated cartography, 
are isarithmic displays shov;ing; for Baltimore tlici distribution of 
such complex phenomena as surface, albedos, energy absorbed and 
emitted by the surface and the distribution of not radiation, often 
considered the injection of energy in a terrestrial tliermal system. 

The project comprises an advance in extraction of data from calibrated 
electro-optical scanners as v.’eil as providing synoptic displays 
of the diversity of energy - related phenomena which can be used 
as analytic tools^ for investigating surface energy states. 

2. Early in January, 1975, the Geography Program received tapes of 
thermal Sky lab data made over the CAI^STS test site on August: 5, 1973. 

The problem was presented to (1) format and interpret the values ■ 
for pixels of the satellite thermal scan and to ( 2 ) correct the 
values recorded for the gray-window affects of the atmosphere and 
thus achieve true surface values. The gray-v;indow in effect damps 
the signal so that a 27®C target records in this case only as about 
2® to 4®C. 

Formatting and deciphering of the hexadecimal tape code was accomp- 
lished by P\obert McNel.ly of the computer center. His ' results were 
mclc according to MASA/BK7LAB directions and were verified by the 
fact that rational energy values were achieved v;hich could be 
calibrated to knovm surface values by the following procedures. 

The gray window ph;sical model, upon which surface to sensor 
calibration was finally accomplished, was originally created as a 
part of USGS/U.C. Riverside participation in the Barbados (BOllEX) 
experiment. In it, a surface signal is Diodified in two ways by air 
through v/hich it is transmitted; (1) the air, by its water vapor 
and aerosol content, alternates or V7eakens the signal, but (2) this 
. loss is to an extent compensated by tha addition of up’ward atmospheric 
longwave emission. Only the attenuated surfcice signals can produce 
a thermal image, the air addition acts as unwantG'd, contrast-reducing 
noise. In its original form, the model was handled as a single equation 

For purposes of calibrating Skylab thermal pixel data, the gra^^- 
vindow model was brnKen into its tv7o parts whi.ch w-ere considered 
separately. This arises from the fact that the satellite is in space 
and the air emission is actually simply atmospheric longwave radiation 
•. leaking into space. On the other han4, the signal attenuation is 



Of course , both 


takinj^ place in the lowest atncspheric layers, 
factors are eventually brought together. 

It was the original intention to estimate the upward lon?y.^ave ^ 
atmospheric radiation (R^a) from the amount of water vapor rn^the 
air chunn, as Masured in precipltable ccntiretcrs as derivea iron 
a radiosonda sounding cade at Baltir=ore Washington Internation Airport. 
The air colunn was to be cj:ar.iined in its stratification using 
follox^ing equation for each layer: 

Py = m (P^ - P2) / S 

Where = Water Vapor in precipitable centimeters, n 

ratio of each layer, (?i - ?■.) Che mean air pressure of the layer in 

dvnes a nillibar= lOOC dynei) , and g the acceleration 0 . ' 

980 cm/sec^. Some nroblen still exists as to how to translate the 
precipitable water value into Ria since values must be for radiance 
in the spectral sensing wlndov; only. Determination 0 . water vapor 
transmissivity by Beers Law (e"i^n) nay be best. 

Since the atmospheric soundings were not yet available, an alternate 
method was devised to determine both precipitable water and atmospheric 
energy emission to space. First, synoptic weather conditions ‘he 

overflight time were determined from appropriate weather maps whic . 
with the im.agcrv defined the depth of the humid laver and its dew- 
point. By the above equation, tiie water vapor in cuis layer was giv 
hose estimate. Second, a numerical iteration was performed in which 
various values for RTa were inserted into appropriate equations (see 
attachment 1) and by means of the Beer's Law relationship (e ) , 

: or ronhnts were'determined to match each v^ue of atmospher e upward 

emission. The value that matched the estimated water vapor content was . 
dLmed to be close to correct. Thus, values derived by two different 
approaches must agree. 

To enable this iteration process to work, a new refineirent was 
added to the gray-window model. Beer’s Law as used app les on y 
v^ater vapor in the air. A method was therefore cevised to measure a 
complementary tuibidity transmissivity (It)* The overall transmissivi y 

(T) thus became: 

T = Tt VThere T„v == water vapor transmissivity 

This constitutes a worthwhile addition to the model which should make 
it better fit reality. . 

It is possible that the numerical solution is as improvement on _ 
original plans for a calibration lacking calibration targets. Since 
the original grmind truth target was omitted, data was obtained tha 
enabled Chesapeake Bay to act as a substitute. This one target made 
it possible to tie calibration equations to the surface. An anaiy. c 
solution for use when one calibration target iS a-zailabU is included 
as attachment 2. The basic c((uatien for the sensor to surface correet 

is : 



Rq = where = spectral radiance at surface 

T 

R^ “= spectral radiance measured aloft, R'^a = 

Long\ 7 ave atmospheric emission as a spectral radiance in the sensing 
window, 

A successful tentative calibration was achieved by the procedures 
described. These permit computer programs to be written to convert 
Skylab thermal tapes into line-printed graymaps showing actual surface 
radiation temperature distributions at the time of imaging. The 
calibrations will be further checked when the atmospheric soundings 
for August 5, 1973 are available from Asheville, A procedure was 
provided to guide the programmer, 

Black/v 7 hite transparencies for the thermal band of the tape plus 
three bands made with reflected solar energy have been received. These 
can be made into a sequence of energy - related maps as v;as done for 
Baltimore if desired. Such a set of maps, with the broad geographic 
perspective of a satellite field of viev.' v;ould indeed be a first and 
should be very useful in the analysis of CARETS land-use climatology. 

The map sequence can also be made directly from the original tapes 
with proper programiming. 

Success of Skylab calibration suggests that satellites are feasible 
platf orms for thermal scanning and provide a much broader geographical 
field of view than is possible with airborne platforms, 

e, Smma r y outlook for the remaining effort to be per formed : No change. 

f. Travel summary a nd plans : Robert Pease visited the US6S liJational 

Center from January 15-17 to analy?;e the tape data as described in 
section d. 
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ATTACH'fENT 1 


Equations Used to Kumerically Iterate for Einission to Space and 
Precipitable VJater 

Kj, = spectral radiance measured at 
satellite in \;/cm‘'/ster/ym 


1» Rz “ Ria 


R 


o 


= T where 


= spectral radiance at surface in 
w/cm^/ster/ym 

Rta = spectral radiance of atmosphere 
into space in w/cm^/ster/ym 

T = overall transmissivity of air 
column 


2. T = T^(\tv) 


T = turbidity transmissivity 
t 


where 

T = water vapor transmissivity 
vrv 



= Beer’s lav? exponential expression 
of v;ater vapor transmissivity 

k = kondrateyev mass absorotion 
coefficient for 10-12 ym 
spectral sensing v:indow = .10 


u - water vapor in precipitable 
centimeters 



Silicon cell measured sol ar input 

Theoretical solar input (with sun angle and water vapor correction) 



ATTACH>OiNT 2 


An Analytical Method for Determining Rfa V^hen One Target is Known 


1. Determine precipitable water in air column by equation: 

= water vapor in precipitable centimeters 
m = mean mixing ratio 

(’’l - P2) = mean pressure of air layers in dynes 

= in(Pi - ?2) 2 

^ g where g = accelleration of gravity = 980 cm/sec 


2. Use P^ = u in Beer’s Law (e ) exponential expression for water 
vapor transmissivity (T^^) , K = mass absorption coefficient 
(kondrateyev)= .10 for 10-12 pm. Resolve the exponential to deter- 
mine the for u amount of precipitable water, 

3. Find overall transmissivity T = x Determine (turbidity 

transmissivity) as described, 

4. Multiply knowii target Rq by T and subtract product from R^ to find 
Rfa. 


R, - R 'I == Ria 
z o 

This becomes a constant for the image. 

5. To solve for other image elements: 


T 



